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Abstract 

A new phase field model is introduced, which can be viewed as a 
nontrivial generalisation of what is known as the Caginalp model. It 
involves in particular nonlinear diffusion terms. By formal asymptotic 
analysis, it is shown that in the sharp interface limit it still yields 
a Stefan-like model with: 1) a generalized Gibbs-Thomson relation 
telling how much the interface temperature differs from the equilib- 
rium temperature when the interface is moving or /and is curved with 
surface tension; 2) a jump condition for the heat flux, which turns out 
to depend on the latent heat and on the velocity of the interface with 
a new, nonlinear term compared to standard models. From the PDE 
analysis point of view, the initial-boundary value problem is proved 
to be locally well-posed in time (for smooth data). 
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1 Introduction 



Phase field models are widely used in various physical contexts in which a 
material exhibits two distinct phases. This is the case for solid-liquid mix- 
tures (e.g. ice-water or alloys during solidification) or for liquid-vapor mix- 
tures (e.g. boihng water), but also for elastic materials subject to martensitic 
transformations. The phase field approach is of special interest, in particular 
for numerical purposes, when interfaces between the two phases are expected 
to show complex geometries and topological changes. In phase field models, 
the 'interfaces' are actually viewed as diffuse interfaces (see for instance the 
famous review paper [T]), i.e. transition regions of nonzero thickness across 
which a so-called order parameter varies smoothly from one to the other of 
its values in the distinguished phases. Here we are interested in a phase 
field model designed for solid-liquid mixtures at rest, which consists of an 
AUen-Cahn type equation for the order parameter coupled with a modified 
heat equation taking into account both the latent heat and the increase of 
entropy due to the non-equilibrium situation inside phase-transition regions. 
This model turns out to be a refined version - in a nontrivial way - of what is 
known as the Caginalp model [3], and it can also be viewed as a special case 
of another one designed by Ruyer [T3j for moving liquid-vapor mixtures. 

The aim of this paper is twofold: 1) by formal asymptotic analysis, we 
show that in the sharp interface limit our model yields a Stefan-like model 
with a (generalized) Gibbs-Thomson relation telling how much the interface 
temperature differs from the equilibrium temperature when the interface is 
moving or/and is curved with surface tension, together with a jump condition 
for the heat flux, which turns out to depend on the latent heat and the 
velocity of the interface with a new, nonlinear term compared to standard 
models; 2) from the PDE analysis point of view, we prove the local well- 
posedness of the Cauchy problem and initial-boundary value problems for 
smooth data. Given that our model displays nonclassical features - it may 
be seen as a degenerate react ion- diffusion system with nonlinear diffusion - 
global well-posedness or rough data are not addressed here. 

The mathematical literature on phase- fields equation is extremely vast. 
In particular, there exist many extensions of the original Caginalp model 
developed in [3] . Let us in particular refer to [lOl |2l |6l [TTl |5l [9l [12] , which are 
quite recent papers, and to references therein. We will not attempt to give 
a minute comparison between these references and our work: let us simply 
emphasize that, up to our knowledge, the model we consider here is distinct 
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from all the models considered so far, mainly by the occurrence of the second 
order quadratic term Ay:>dt(p {ip being the order parameter) in the equation 
for the temperature, cf. the second equation in f|T6l) . 

The paper is organized as follows. In Section [2] we derive the model and 
its six-parameters nondimensionalized version. The sharp interface limit is 
investigated in Section [SJ Local well-posedness is shown in Section HJ 

2 Phase field equations 

2.1 Derivation and basic properties 

The model we are going to consider pertains to the so-called second gradient 
theory. We assume that the physical state of a solid-liquid mixture is de- 
scribed by an order parameter ip and its temperature T in such a way that 
its free specific energy / depends on T, ip and also Vip in the following way 

(1) f{T,p,Vp) = ^ (w{p) + ^\\Vp\'^ - j\siT,p)dT, 

where p is the density of the mixture, which will be assumed to be homo- 
geneous and constant, Tg is the equilibrium temperature, A is a positive pa- 
rameter that is supposed to govern the width of solidification/melting fronts, 
is a double-well potential, and s is the specific entropy of the mixture. 
More specifically, the order parameter is chosen so that in the pure phases 
we have either p = 1 (liquid) or = (solid), and W is supposed to achieve 
its global minimum at both and 1 and nowhere else. Furthermore, s is 
taken to be a convex combination of the entropy in the phases, depending 
nonlinearly on the order parameter in the following way 

(2) s{T, p>) = u{p>)sn^{T) + (1 - i^{p))s,oi{T) , 

where u : [0, 1] — )■ [0, 1] is monotonically increasing. Typical graphs of the 
functions W and i^' are represented on Fig. [H By contrast, in the phase field 
model of Caginalp [3], u would be the identity function (hence z/' = 1). 
The latent heat of the phase change is by definition 

(3) C{T) := T(sHq(T)-Ssoi(T)). 
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Figure 1: Derivatives of the double- well potential W and of the entropy 
coefficient u. 



So another way of writing the free energy is 
(4) /(T,y,,Vy.) 

and thus the (standard) chemical 'potential of the mixture is 



(5) 



/i(T,(^) : = 



^ C{r) 



dr. 



r 



Since W has wells at and 1, we see that yU = in both phases whatever the 
temperature, provided that u' vanishes at and 1 (as on Fig. [T]): this would 
obviously not be the case for a linear u, as in standard phase field models. 
The heat capacity of the mixture is 

(6) C, := T— = z/(¥.)C,,Hq(T) + (1 - u{v))C,MT)- 

To simplify the analysis, we shall assume that the heat capacities of the 
liquid Cp^iiq = Tdsn^/dT and of the solid Cp^soi = Tdssoi/dT have the same 
constant value Co, so that Cp is also equal to Cq. In other words, we shall 
concentrate on the special case 



(7) Cp = Co, C{T) = C,- 
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Regardless of that simplifying assumption, we consider the following equa- 
tions for the evolution of the mixture: 



(8) 



pCpdtT + p [pg-T^^ dt^ = div(WT) 



where /c > denotes the heat conductivity, k > denotes the so-called 
mobility, and 'generalized chemical potential', which merely differs 

from the standard chemical potential by a second order term: 

(9) /x,[T, V'] := ^ = V.) - ■ (AV^) . 

Observe in particular that /i^ = in the phases = or = 1), as for 
p,. The first equation in ([HD is the building block of phase field models, in 
which 1/k is presumably proportional to a relaxation time for the mixture 
to return to equilibrium. Taking into account that 

^ dT' dT dT dif' 
the equations in (|H]) ensure that (for smooth solutions) 

(10) pdts = div -—- + k ' ' + ^ 



V T / T 
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which (formally) means that the growth of total entropy / s{t,x) dx is gov- 
erned by both the conductivity (k) and the mobility (k). 

In fact, dH]) is specifically designed to have (fTOj) as well as the (formal) 
conservation of total energy. More precisely, the specific energy e = f + sT 
is conserved along solutions of ([S]) in any domain f2 C M'' such that 

(11) / XVif-ndt^ = and / kVT ■ n = , 

Jan Jon 

where n denotes the normal to dQ. Indeed, recalling that s = —df/dT, the 
first equation in fITT]) enables us to write 

— [ edx = [ ^dt'fdx + [ Tdtsdx, 
dt Jn Jn o^p Jn 
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where, by definition of fig and by tlie first equation in dH]), the first integral 
equals — f hi fi"^ dx, which obviously cancels out with the integral coming 
from the last term in Eq. ffTOj) . To conclude that e{t, x) dx is constant, we 
observe that by the condition on T in ffTTl) . 

To finish with these general observations, we point out that the equali- 
ties in ( ITTl) are easily achieved by means of standard boundary conditions. 
Namely, the second equality will be implied by a homogeneous Neumann 
boundary condition on T, i.e. VT ■ n = (meaning zero heat flux at the 
boundary: incidentally, one may note that for a nonzero heat flux the total 
energy will either decrease, due to cooling, or increase, due to heating), and 
either a homogeneous Neumann boundary condition or a Dirichlet condition 
ip\g^ = or 1 (both values implying fig = 0, as already noticed) will ensure 
the first one. The appropriate choice of a boundary condition for ip is related 
to the moving contact line problem, which we shall not discuss here. 

2.2 Nondimensionalization 

The total number of independent physical units used to describe all depen- 
dent variables and independent variables in ([8])-([9]) is ten (those of x, t, k, 
W, p, T, £, A, Cp, k, and of course ip and u do not count because they are 
already nondimensional) , and the number of fundamental physical units is 
four (kg, m, s, K). So by elementary dimensional analysis (Buckingham tt 
theorem), a nondimensionalized version of ([5])-([8])-([9]) requires 10 — 4 = 6 
nondimensional parameters. Below is a possible choice for these parameters, 
expressed in terms of 

• the density p, 

• the equilibrium temperature Tg together with a characteristic temper- 
ature difference 6T, 

• a length scale L, 

• a characteristic interface thickness h, 

• a time scale Iq, 
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• the surface tension a, 



• the latent heat £e at Tg, 

• a reference heat capacity Co, 

• a reference mobihty coefficient kq, 

• a reference heat conductivity /cq- 
Introducing the parameters 

e ■= - Pe •= P^^^^ a ■= 9 ■=— P ■= ^ St •= ^^^^ 

L ' koto ' ph ' ST ' pCoh6T ' 

together with the rescaled variables 

~ 2; ~ t ^ ~ W ~ T-Te ~ £ ~ A 

x:=-, t:=-, K := — , := — - , T := , £ := — , A := — , 

L to kq a/n 01 Le crh 

-Co'" - ko ' 

we may rewrite (IE])-© as 

d-jip = — aHJlg , 

^^^^ ^ c,d~r + p (^Jl, -{9 + T)^) = i- ■ (W,f ) , 



(13) J1,[T, ^] = W'{^) - ^ u'{^) ^ ^ dr - 5^ V, ■ (AV,^) , 

If K is supposed to be constant, or similarly if k is constant, we may assume 
without loss of generality that k = 1, respectively /c = 1, in (fT2|) . The case 
of A in f lT^ is more subtle because it depends on the parameters h and A, 
the former being arbitrary and the latter not being accessible to physical 
measurements. Nevertheless, we choose to set A = 1. In addition, under the 
simplifying assumptions in ([7]), we have in the rescaled variables 

= 1 , £ = 1 + ^ , 
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so that in this case the nondimensionahzed version fll2l) - fll3l) of ([H])-® reads, 
dropping the tildes for simphcity, 

dt^ = -a fig, 

^^^^ ^ dtT + P (^fi,-{e + T)^'^ = 1-AT, 

(15) /x,[T, ^] = W'icp) - ^ i^'icp) j-e'Acp. 



Plugging (1151) into flMj) we get the system 
(16) 



a T 



This resembles the system considered by Caginalp in his seminal paper [3], 
except for two important differences. The first one is that the coefficient of 
T depends on (f in the first equation. The other one lies in the complicated, 
second order and nonlinear coefficient of dt^p in the second equation, which 
is - up to the authors knowledge -, always supposed to be a constant (latent 
heat) in Caginalp-like models. 

For completeness, let us now derive the nondimensional versions of the 
entropy equation f lTOj) and of the local conservation law for the energy. Re- 
defining s as the nondimensional entropy s/Co, we have from that 

' = ^^(^) + HT + e) 

up to a harmless additive constant. Then the nondimensionahzed version of 
the entropy equation (fTOj) is 



Pe \T + eJ Pe {T + ef a{T 
Regarding the nondimensionahzed energy 

e:=/+ {T + e)s = T + (^W{^) + ^ z/(v^) + ^/3e'|V^P 
we easily find the conservation law 

(18) dte = ^AT + /3e^V- {{dtif) Vif) . 

Pe 
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3 Sharp interface limit 



Our aim here to derive at least formally a physically realistic, asymptotic 
limit of the system ( !T6|l when the width of interfaces tends to zero, either 
because of a physical scaling or for other reasons related to the actual values 
of the six nondimensional parameters a, (3, e, 6, Pe, and St. More precisely, 
we are going to show in what follows that for suitable relationships between 
those parameters, the system f|T6|l formally tends to the Stefan-like model 
f l3T|) (see p. [17] hereafter) when e goes to zero. (Recall that e is the parameter 
governing the typical width of interfaces.) Before entering into details, let us 
emphasize that the sharp interface model in f lHT]) naturally involves the heat 
equation in the phases, and two sorts of conditions at interfaces, namely 

• a (generalized) Gibbs-Thomson relation giving the interface tempera- 
ture in terms of the surface tension, the mean curvature and the velocity 
of the interface, 

• a jump condition for the heat flux across the interface, in terms of the 
latent heat and of the velocity of the interface, the later dependence 
being nonlinear (quadratic). 

This should be of interest to discuss the physical validity of f ll6p . 

3.1 Formal asymptotics 

For convenience, we rewrite f ll6p as 

adtv = e^Aif - W'{ip) + ^v\^)T, 
PdtT = 6 AT - ^ {T + 6) dtu{ip) + a{dtip)\ 



(19) 
with 



1-^1 1 1 

(20) «:=-, ^■■=^. l--=l^n^ S:- 



' ' (3St9' /3Fe 

The six nondimensional parameters in (fT9l) are now a, 13, 7, S, e, and 6. If 
we go back to the original definitions of a, /3, St, 6, and Pe, we see from (l20!l 
that a, /3, 7, 6 are all proportional to the ratio h/a of the interface width 
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and the surface tension, and each of them has its own a physical significance 
according to the following relationships 



d oc I/kq (relaxation time) , 

(3 (xCq (heat capacity) , 

7 oc £e (latent heat) , 

5 (xkQ (thermal conductivity) . 

As regards the sharp interface limit e = h/ L ^ Q aX fixed surface tension a, 
by the observation above it is rather natural to let the four parameters d, /3, 
7, and 5 go to zero at least like e. If in addition we let the relaxation time 
go to zero like £, we are led to consider 

a := a/e^ , /3 := P/e , ^ := 'j/e , 6 := 5/e 

as being fixed. With these definitions, fll9p becomes 

ae^dt^ = e^Aif - W'{ip) + ^ev\^)T, 



(21) 



The formal limit of the first equation in ( l2T]) as e — >■ gives W'{(p) = 0, 
which imposes that ip takes only the values (solid phase), 1 (liquid phase), 
or a ('metastable' state), while the formal limit of the second equation is 

(22) ^dtT = 5AT - -f {T + e) dMip) . 

Assume that T is a continuous solution of (l22l) . in which ip represents a 
sharp interface, that is, ip is constant and equal to or 1 on either side a 
smooth, moving surface r{t) as on Fig. [2l Then by integration by parts in 
the neighborhood of any point (x, t) G S := {(x, t) G M'^ x M , a; G r(t)} we 
find that the gradient of T experiences a discontinuity across r{t) according 
to the following relation 

(23) -^{T + e)v = 6 [VT ■ N] 

where denotes the unit normal to F pointing to the liquid phase (</) = 1) 
and V denotes the speed of F in the direction A^. A linear relation between 
V and [VT ■ A^] as in Equation (l23l) is a classical building block in Stefan 
models for sharp interfaces, see for instance Fig. 1 in [4]. 
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Figure 2: Sharp interface configuration 

Of course the formal hmits above are not vahd in regions where (y?, T) 
experience large variations, and to describe exact solutions of (12T]) we need 
internal layers for diffuse interfaces. In what follows we adopt the same, 
multiscale approach as in [1], where the sharp interface limit was obtained for 
the usual Caginalp model. Consider a (smooth) solution (c/^^, T"^) of and 
let T'^it) be the level surface {x G M*^ ; ip'^{x, t) = b} (the b where u' attains its 
maximum, which is supposed to best describe the location of the 'interface'). 
We assume that T'^{t) is smooth, not self- intersecting, and depends smoothly 
on t and e in such a way that the signed distance d'^{x,t) of x to T^(t) is 
well-defined for t E [0,T], e E [0,eo], and x in some neighborhood 'fit) of 
UeG[o.eo]r^(t). By definition of d'^, 

N'{x,t) := Vd'{x,t) 

is a unit normal vector to r^(t), 

v^{x,t) := —dtd'^{x,t) 
is the normal speed of r^(t) in the direction A^^, and 

H%x,t) := -Ad'{x,t) 
is the sum of principal curvatures of r^(t) at x. Assume that the functions 
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d^, (f^, and T'^ admit asymptotic expansions of the form 

oo 

t) ~ ^ 5* di{x, t) , 

oo 

i=0 

oo 

as e — 0, for x G 'J^(t). We recall that the notation ~ Xli^o^**^* means 
that for all / G N, rf= - E[=o^*^i = We shall denote hj z = 

d^{x,t)/e the rescaled variable in the normal direction to r^(t). For any 
smooth function {x^t^z) H- F{x,t,z), the derivatives of the function : 
{x,t) F{x,t,d'^{x,t)/e) are given by 

dtF' = dtF + e-' {d,F) dtd' , 

VF" = VF + e-^ {d,F) Vd' , 

AF" = AF + e-^ {{d,F) Ad' + 2Vd' ■ V9,f) + (9,^^ , 

where the differential operators V and A concern only the variable x. Hence 
the system fl?Il) for {^p,T) = ((y9^,T^) is equivalent to the following one, 
evaluated at z = d'^{x,t)/e: 
(24) 

' (9,V - W{if) = e ({adtd' - M')d,^ - 2Vd' ■ Vd,^ - 7 z^'(^) 
+ {adt^ - A^) , 
Sd^f = e {(pOtd' - dAd')d,f - 25V d' ■ Vd,f 

+ ^{T + ey{^){dtd')d,^ - 2a{dtd'f{d,^f) 

+ {^dtf-~5f + ^{f + e)dtu{^) - a{dtd%dt^)d,^^ 

We expect that r^(t) converges to r'^(t), the level set {x G M*^ ; do{x, t) = 0}. 
Off r°(t), d'^{-,t)/e tends to ±oo as e -> 0, so that we shall need extensions of 
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T) for all values of 2; G (—00, +00). However, the only constraint is that 
f l24p holds at z = d'^{x,t)/e, which means that we can add to the equations 
any 'reasonable' function of {x,t,d'^{x,t) — ez) that vanishes when its last 
variable equals zero. This observation will be used in a crucial way to deal 
with the equation on T. 

Retaining only the £° terms in the first equation of f l2^ we get 



(25) 9,Vo - W'{^,) = 



which is the standard equation for a stationary diffuse interface connecting 
at 2 = —00 to 1 at z = +00 (or vice- versa). A straightforward phase 
portrait analysis shows that there is a unique such tp^ satisfying v5o(0) = h. 
In particular, there is no degree of freedom for tpQ to depend on (x,t). 

To the next order, using that (p^ is independent of x, we obtain from the 
factors of the equation 

(26) {dl - W"{^o)) ^1 = {adtdo - Ado) d.^o - wivo) ^ . 

Since ipo tends to its endstates exponentially fast, the right-hand side of 
fl26|) tends to zero exponentially fast provided that To is bounded, or has at 
most polynomial growth in z. In this case, since by differentiation of ( l25l) the 
derivative of the interface profile dzfo belongs to the kernel of the self-adjoint 
operator 9^ - W"{(po) in L^{R) (with domain H'^{R)), a necessary condition 
for (!26l) to have a solution ipi{x,t, ■) G H'^{M.) is 

/ + OO P + OO 

{d^ipofdz = 7 / TQd,u{ifQ)dz . 
-00 J —00 

Defining the 'interface temperature' by 

(To) := / ndMVo)dz = ±^"~ 
J —00 



-00 , 
-00 
"+00 
-00 



/_r dzi^i^o) dz 



and the (nondimensional) surface tension by uo := J_^{dz'PoY ^^e previ- 
ous relation may be seen as a (nondimensional) generalized Gibbs-Thomson 
condition: 

(27) ao (adtdo - Ado) = ±1 (To) ■ 

Here above the ± sign is merely a shorthand for [^^(v^o)]^^, which equals 
if N^{x,t) := Vdo{x,t) points to the liquid phase (or —1 of points to 
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the solid phase). RecaUing that Vq := —dtdo is the normal speed of and 
-f^o ■= —^do is the sum of principal curvatures of r°, we can indeed identify 
f l27|l with the usual condition in generalized Stefan models (see again Fig. 1 
in i). 

As regards the second equation in (^^, the only term of order zero in e 
is 6 d^T. Nevertheless, we may add to that equation a function of the form 

h^{x, t) p{z) (c/^(x, t) — ez) , 

which obviously vanishes aX z = d'^{x,t)/e, with p smooth and compactly 
supported in z (so that the term ez is at most of the order of e) and h^{x, t) ~ 
/ij(x, t). More precisely, we shall assume, similarly as in |3], that 
p = dirj with Tj such that 77 = on (—00, —1], 77 = 1 on [1, +cxd), and 77' > 
on (—1, 1). Then the zeroth order equation becomes 



5diT^ = hoix,t)do{x,t)dir], 

which necessarily yields, if To is sought bounded in z, that 

6To{x,t,z) - ho{x,t)do{x,t)r]{z) =: 6To~{x,t), 

a function of (x,t) alone. Since ri{z) = for z < —1, this is a consistent 
notation in that TQ{x,t) = \imz^^aoTo{x,t, z). Moreover, since r]{z) = 1 for 
2; > 1, we have 

ho{x,t) do{x,t) = 6(T^{x,t) - To"(x,t)), 

where Tf^{x,t) := lim^^+oo 7o(x, t, 2). This shows in particular that for ho 
to be smooth, Tq and Tq" must coincide on the zero level set of do, namely 
on r*^. Conversely, if To*" and Tq" are smooth functions coinciding on r°, we 
can define 



ho{x,t) :-- 



^ To+(x,t) T^{x,t) ^ x^ro(t) 
do{x, t) 

5 [VTo ■ iV°] (x, t) , xeT^it) 



where = Vdo (as before) and the 'jump' notation [VTo ■ A^°] merely stands 



for VTo+ ■ iVO - VTq- ■ A^^. Then, 



To{x,t,z) := 6 ho{x,t) do{x,t)'r]{z) + Tq (x,t) 
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is independent of z for x G r°(t), and more precisely, 

To(a;,t,2;) = To~(a;,t) = To+(x,t) = ±(To)(x,t), 

where again ± = [i^(v5o)]i^- 

Now, the next order terms in the asymptotic expansion of the second 
equation in fl2^ supplemented with the term h'^{x,t) {d'^{x,t) — ez) d'^rj^z) 
will enable us to find a necessary relation between the heat flux [VTq ■ A^*^] , 
the interface temperature (Tq) + 6, and the velocity = —dtdo of T^. As a 
matter of fact, retaining only the terms of order one, we get 

Sd'^Ti = {hido +hodi)d'^r] — hozd'^r] 

+ [Mdo - SAdo)d,To - 26Vdo ■ Vd^T, 

+ 7(To + ^) {dtdo)d,u{^^) - 2a{dtdof{d,^Qf. 

A necessary condition for Ti to be bounded in z is that the integral from z = 
— oo to 2; = +00 of the right-hand side above equals zero. The contribution 
of the first row to the integral is just /iq, because / dlrj = and / zd^r] = — 1 
by definition of rj. The next term does not contribute if we restrict to a; G 
r°(t) because then To+(x,t) = TQ{x,t). Recalling that ho = 5 [VTq ■ N°], 
Tq = ±(To) on r°, and that we have defined 

/+00 
-00 

we finally arrive at 

(29) 6 [VTo ■ N'] = ±7 ((To) + e) (dtdo) - 2aao {dtdof . 

Observe that the roughly obtained relation fl23|) may be seen as an approxi- 
mation of ( l29i) when the velocity = —dtdo of r° is small enough. 

To summarize, the sharp interface limit of (l2Ti) is expected to be (rigor- 
ous justification will be addressed elsewhere) the generalized Stefan problem 
consisting of the heat equation for T outside r° together with the conditions 
(!27|) . (!29l) on r°. If for instance A^ = points to the liquid phase, this 
problem reads 

( pdtT = 6AT outside r, 

(30) < aoiH-av)=^T on T , 
[ 6 [VT ■ N] = -^{T + e)v - 2aaov^ on T , 
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with (Jo defined in fl28l) where ip^ is solution of (123]) and tends to at — cxd 
and 1 at +oo, v the normal velocity of F = r°, and H the sum of principal 
curvatures of F. 



3.2 Back to physical variables 

The sharp interface model (15U]) is non-dimensional. It is of course important 
from the physical point of view to go back to physical quantities. 

Let us start with the first equation in ( 130|) . which by definition of (5 and 
5 also reads 

Y>edtT = VT. 
Remembering that T, t, and x respectively stand for 
~ T — Tf, ~ t 

(where the tilda are those of Section 12.21 and not those of Section 13. ip by 
definition of the Peclet number Pe = pCpL"^ / {kto) we recover the expected 
heat equation 

pCpdtT = kAT. 
As to the last equation in ( l30l) . it actually reads 



1 

Pe 



— -(T + 0)v — 2 V 



with V = &jd/L = {tQ/L)v if v denotes the physical velocity of the interface. 
Before going further, let us comment on (Tq = f^^idzfoY dz, where ipQ is by 

definition (see Eq. fl25l) . having in mind that W stands for W = {h/a) W) 
solution of the differential equation 

An obvious integrating factor is d^ipo, and since W{ipo) vanishes at ±oo, we 
have {dz^Po)"^ = 2{h/a) Wi^ipo). This implies that 

2 r+oo 

ao = - W{ipo)hdz. 

^ J~oo 

Recalling the meaning of the z variable, which scales as the actual distance to 
the interface over eL = h, we can identify the integral 2 W{(po) h dz with 
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the physical surface tension a, and therefore set ctq = 1. Substituting all the 
other non-dimensional parameters Pe, St, 6, /3, a and e by their expressions 
in terms of physical quantities, we get in turn the physical jump condition 



[VT ■ N] 



Cv 



K, h 



(Recall from ([7]) that C = CeT/T^.) This is to be compared with the usual 
jump condition in Stefan models: 



P 



[VT ■ A^] 



In particular, we observe that the quadratic correction in the velocity v is 
neghgible if the velocity v is small compared to Cnh (which is indeed homo- 
geneous to a velocity). 

We finish with the derivation of the generalized Gibbs- Thomson relation. 
Recalling that we have set do = 1, that v stands for v = {tQ/L)v and noting 
that H stands for H = A^d = LAd, we can rewrite the second equation in 
([30]) as 

a-^f — LH = — 7T, 
Ij 

where H := Ad is the actual sum of principal curvatures (homogeneous to 
the inverse of a distance). Substituting a and 7 for their expressions, this 
eventually gives 

P 



Kh If 



P^{T 

e 



Therefore, in physical variables the (generalized Stefan) sharp interface 
model ( 130|) reads 



(31) 



{ pCpdtT = kAT 

- [VT-iV] = - 2-- 



outside r 
o H on r , 
on r . 
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4 Well-posedness 



We now turn to the mathematical analysis of the (non-standard) PDEs sys- 
tem (fT9l) . which we equivalently rewrite as 




a dt^ - e^A<^ + W'{<^) = 7 T, 

^dtT + ^e u\ip) dtip -6AT = F{ip, Aif, T), 



with 

(33) F{ip, Aip, T) := \{e^A^ - W\^)f + l{e''A^ - W\^))u\^)T. 

a a 

In this system, d, /3, 7, 5, e, and 6 are fixed, positive parameters, and the 
functions W , v are supposed to be nonnegative and to belong to (M) (the 
space of functions that are bounded as well as their derivatives up to 
order 3). We are going to show that the Initial Boundary Value Problem for 
f p2|) with suitable initial and boundary data is locally well-posed both in two 
and three space dimensions. 



4.1 Functional framework and main results 



In what follows, is an open, bounded, and regular subset of M*^, d G {2, 3}. 
As boundary conditions on we consider a homogeneous Neumann con- 
dition for the order parameter y?, and a mixed constant Neumann-Dirichlet 
boundary condition for the temperature: 



(34) 



on 9^2, 



gb on r, T = Tb on dVt \ F, 



where F is a given, relatively open subset of c?f2, and n denotes the unit out- 
ward normal to c?f2. We suppose that gb and Tb are constants, corresponding 
respectively to the heat flux and to the temperature imposed on the bound- 
ary of the domain. Given gb? ^b? we know from [SI Notes of chapter 8] that 
there exists T e if^(fi) solution of 



(35) 



f AT = in 



On 



Qh on F, 



^ T = Tb on 9fi \ F, 
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in the sense that T - Tb G H^{n U T), the closure of <4^(fi U T) in H^{n), 
and 



VT-Vrdx = J q^rd-W^-U 



x] 



for all r G -f^o(^^ U F). Two solutions Ti and T2 of the first two equations 
in (153]) differ by a constant, so that T is unique if \dQ \T\ > 0, and unique 
up to constant in the case of a pure Neumann condition. 
For s > 0, we denote by H^{Q) the closure of 



..Pin); 1^ 



= 

an 



in H^{Q). In particular, H^{^) is merely the set of functions (f G H'^{Q) 
that satisfy the homogeneous Neumann boundary condition dn^p = on dfi. 
When no confusion can occur we shall just write for for H^{Q), 

and for //^(r^UF). 

Definition 4.1 (Weak solution) Fory^o G if^, Tq G i^o^ ^ (O^+oo]; 

we say that {(fijT) is a weak solution of ( l32l) on [0,t*) wi/i initial data (fo, 
Tq, and boundary conditions (1341) . z/ 

y.G'^,([o,r);i/2)nLL(o,r;i/,^), 

T - f G ^,([0, r); n lL(o, r; i/i), 

wt/i (/5 and T satisfying v^j^^g = v^o, ^l^^o ~ -^O' ^'^^ 

• i/ie /irsi equation in ( l32i) m i/ie sense i/iai 

(36) « /7 ¥'5*C= /7 /7 z/'MTC-e' /7 AcpC 

Jo Jn Jo Jn Jo Jn Jo Jn 

for allCe^o\{OXy,L^), 

• the second equation in ( !32l) in the sense that 

(37) 

f [ 0T + -feu{^))dtT= 6 f [ VT-Vr - / / F(y.,Ay.,T)r 
Jo Jn Jo Jn Jo Jn 

forallr G ^oH(0, t*); n^((0, T); i/g^), where F is defined as m ([331). 
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Note that, according to this definition and by the Sobolev embeddings 
and (both vahd in space dimension c? < 3), a weak 

solution {^p,T) is such that G ^fe([0,t*) x Q), and Aip G L^{Q) at almost 
all times in [0,t*). This gives sense in particular to the last integral in f l37|) 
if we also note that r G L6(fi) at all times in {0,t*) if r G ^((0, r); ifg^): 
indeed, examining all the terms in the product F{ip, Aip, T) t we see, using 
that W and z/' are bounded, that we 'only' need r, Tr, AipTr, Aipr, and 
(A</))^ r being integrable, which is certainly the case on a bounded domain 
when T e L'^,T e L^, and A(f G L^. 

Theorem 4.2 (Existence of weak solutions) For all </?o ^ ^ ? 

there exist t* > and a weak solution {ip,T) of (!32|) on [0,t*) with initial 
data ifQ, Tq, and boundary conditions (!34l) . in the sense of Definition 4-i 



Theorem 4.3 (Continuous dependence on the data) 1. Given (fo ^ 
H^, Tq ^ Hq, there exists at most one weak solution {ip, T) to (!32l) - (!34ll 
with initial datum {ipq^Tq). 

2. If {ipi,Ti), i G {1,2}, are two weak solutions to f l32p -( 15^ both defined 
on [0,t*], and with initial data {(Pifl,Ti^Q), i = 1,2, respectively, then 
there exists a constant C > depending only on the norms 



||^i|U°°(0,i*;//l(^))' ll^i|U2(0,t*;L2(C))) 
||'^i||L°°(0,t*;H2(n)), ||V'i||L2(0,t*;/i'3(n)), 

such that, forO<t<t*, 



^G{1,2} 



(38) ||Ti(t) -T2(t)||i.(^) + Wviit) - Mmhin) 

< C (j\Tifl - r2,o||i2(j^) + ||V51,0 - V52,o||l/2(f^)j . 

In particular, the weak solution to (132|) - (!34|) depends continuously on 
the data. 

Unsurprisingly, the proof of Theorems 14.21 and 14.31 relies on a priori 
estimates. It is to be noted though that we shall use other quantities than 
the total energy 

E{T,^) := [ 0T + j9u{^) + W{v) + le'\V^\^). 
Jn 
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Indeed, if we do have the conservation of E, thanks to fllSp (where e = (3E), 
at least when T satisfies a homogeneous Neumann condition on the whole 
boundary dQ (that is, for qh = and F = dQ), this is obviously not enough 
to control the norm of T. Rather, we shall use 

(39) Eo{T,^) := f {^T' + W{^) + le'\V^\'). 

Compared to -E, the interest of E^ is that it is quadratic in T, and it satisfies 
the identity 

^^Eo(r,v.) + 5||vr||i. + ae\\dt^f^, = j f{^,a^,t)t 

along solutions of (132|) (and (l34l) with gb = 0, Tb = 0). Of course, because of 
the right-hand side, this is not fully satisfactory and we shall need another 
quantity to control the norm of A(y9. 

Remark 1 //we replace F by zero in (152]) . we are left with the Caginalp-like 
model 

\^dtT + ^eu'{ip)dtip-6AT = 0, 
for which we have the much nicer identity 

e^^Eo{T,ip) + S\\VT\\l, + ae\\dM\l, = 0. 

Remark 2 By an adaptation of our a priori estimates, we can show in ad- 
dition that, among the stationary solutions to (132!) . those corresponding to 
single-phase states, i.e. with ip = or ip = 1 and T = T, are stable. 

Remark 3 Since we are interested in asymptotic models, we have chosen on 
purpose to keep track of the nondimensionalized numbers a, (3, 7, 5, e, and 
9 in our a priori estimates. We shall also pay attention to the occurrence of 
\\T\\Hi, and of the bounds for W , v, and their derivatives. 
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Further notations. All constants depending only on the dimension d and 
on Q will be considered harmless, and we shall denote by 

any inequality Ai < CA2 where C is a constant (depending only on d and 
Q). As already mentioned, W and u are supposed to belong to "^jf (M). For 
simplicity, we introduce the notations 

i/^ := sup|z/'|, W^:=sup\W'\, 
and similarly for their second and third order derivatives. 



4.2 Existence of solutions 

In this section, we prove Theorem 14.21 by means of a Galerkin approximation. 

Let {</3j}igN* be a set of eigen-functions of the Laplacian operator —A 
with Neumann boundary condition, {v^ijigN* being a complete orthonormal 
system in H^. Let {TjjjgN* be a set of eigen-functions of —A in Q with the 
boundary conditions 

Of. _ 
— ^ = on r, Ti = on dQ\T, 
on 

with {TjjjgN. being a complete orthonormal set in Hq. We seek approximate 
solutions of (Ea - (El of the form 



^"(t) = ^a,(t)^„ T"(t):=T'^(t)+T, T'^(t) = ^ 6,(t)r, 

1=1 i=l 

where Oj and 6j are 'if^, real- valued functions. Defining 

Vn = Span{v?i, ■■■ ,Lpn}, Zn = Span{Ti, ■ ■ ■ , T„}, 
we require that for all C ^ V„, 



(41) 



Jn Jn Jn Jn 
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and for all t ^ 

(42) I pdtTT = -5 [ wT-Wt-^iO [ i^'((^")at(p"r + / F"r, 
Jn Jn Jn Jn 

where 

together with the initial conditions 

^"(0)=Pv„(^o), T{0)=P^,XTo-f), 

where for any subspace Y, Vy denotes the orthogonal projection onto Y. 

Denoting by a and b the vector-valued functions of components and 
bi respectively, taking ^ = y^j in ( BTj) and r = Tj in (142|) for i = 1, . . . , n, we 
obtain ordinary differential equations of the form 

(43) l| = $(a,b), ^ = r(a,b,$(a,b)), 

with $ and T at least on M^" and M'^" respectively (since W and u are 
'^^). Therefore, the Cauchy-Lipschitz theorem ensures the existence and 
uniqueness on some maximal time interval [0, t„), t„ G (0, +oo], of a solution 
(a, b) with prescribed initial data. 

The next step is to derive some estimates in order to show that the 
times sequence (t„) is bounded from below by a positive time and that the 
sequences (v^"), (T ) are bounded in the appropriate functional spaces. To 
simplify the notations we will drop the superscript n. 



Energy estimate. Recalling the definition of Eq in ( l39|) . and taking a 
combination of (1411) with ( = ipn and ( 142|) with t = Tn, we get the identity 

(44) e^Eo{T,^) + 6\\VT\\l, + a 6 WdtifWh = [ FT + ^6 [ u'{^)dt^f. 

'^'t Jn Jn 

By the elementary inequality 

(45) a6 < Aa^ + — , Va, 6 G M , VA > , 
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we have 

7 / i^\v)dt^f <'^\\dM\l^ + ^^WJ'\\T\l.. 
Jn 4 a 

Recalling that the notation a actually means and introducing the new 

simplifying notation 

(46) ^■.= ^^\u'J\ 
we thus obtain from (j44j) that 

(47) ^i?o(T, ^) + 69 \\VT\\l, + |d \\dM\l^ < 9 j^FT + a \\f\\l,, 

where we have used the notation 9 for 1/9. We shall use the same convention 
repeatedly for other positive parameters in what follows. 

The integral of FT in ( l47l) involves trilinear terms in Aip and T. We can 
get an estimate on them by using the following, higher order estimate. 

Lemma 4.4 (second order estimate) Let us define the modified energy 

E,{T,y^) = Eo{T,y^) + le^6afi9 [ \A^\\ 

Jn 

and (to enforce estimates of powers of Ei) 

E*{T,ip) = max{l,E,{T,ip)). 
Let the constants , B^, C^, be given by 

A° := 6ujil9e^, D° := max(l, 9^) 9/3 (1 + fia96) l/je^I + L/^i^), 

B':=P9\\f\\l,, C°:= (/xa + 5^)||f 

where 

a;:=|W^^P, /i := 7^ 



oo I 



u" 

oo 



Then the first n modes (T = T + T,ip) = (T" = T" + T, </?"■) solutions 
to fHT|) - fH2|) satisfy the refined energy estimate 

(48) ^Ei(T, + 5^11 VTlli. + a\\dM\h + e'6fi9 || VA^||i. 

-9 [ F{^, A^,T)T < A'E,{T,^) + D%B^ + El{T,^)f + C', 
Jn 

for all t G (0,t„). 
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Proof: Let us first write 

Ei{T,ip) = Eo{T,ip) + za I \Aip\' 



n 



with the constant z to be determined in the course of the proof. We apply ( l4Til 
with A^(p as a test function. Since we use a Galerkin approximation buih on 
the eigen-f unctions of the operator —A with Neumann boundary condition, 
there is no boundary term in the next integrations by parts. We obtain 

(49) laj\\A^\\h+s'\\VA^\\l = - [ W'i^)A'^ + ^ j u'{^)TA'^. 

^ V ' ^ V ' 

Jl J2 

Using an integration by parts, we get an estimate on Ji because 

c-2 lU/W 12 

jl|VAv.||i.<^||Vv.||i. 



thanks to the inequahty in fHSl) . Another integration by parts gives 

J2 = --f u{ip)VT ■ VAif - 7 / u"{ip)T Vif ■ VAif. 
Jn Jn 

By the Cauchy-Schwarz inequahty, we thus have 

IJ2I < T^^^ollVTlU^llVAyplU^ +7z/:^||T|U2||V(^|Uc.||VA(^|U2, 
hence by using again. 



\J,\ - V llVA(/.||i. < 



1^0 



2 

2 



VT||i.+7z/;^||r|U2||Vv^||Loo||VA(/.|U2. 



e 

To control the L°° norm of Vy^, we can apply to m := V(f Agmon's inequality 
(50) ^ \\u\\h2 , 

which is valid in dimension d < 3 - in the two dimensional case, it just follows 
from the Sobolev embedding H^^'^ L°° and the interpolation between 
and H^, showing that H^^"^ = [H^,H\/2. Using that \\u\\h2 < \\u\\l2 + 
\\Au\\l2, this gives 

IIV^IU^ < ||V^||^/?(||V^||i^/V IIVA^II^/^), 
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from which we deduce that 

l|r|u.||V(/.|U^||VA^IU. < ||t|u.(||Vv^||hi||va^iu. + ||Vv^||^/,'||vA(/.||^/') 



Then, using once more ( H5l) together with its more general version, Young's 
inequahty 

aP , b'^ 11 

(51) a6< A— + A"«/P-, Va,6eM,VA>0, Vp, g > , - + - = 1, 
p q P Q 

with p = 4/3, we obtain 

\J2\ - J llVAv^lli. < /ielVT||i.+./ie2||T||i.||V^|||,+6V^lT||i.||V^||^„ 
where we have used the announced shorthands 



2 I / |2 



oo 



Now, since > by assumption, we have 

WV^Wl^ < 2e'Eo{T,^), \\T\\l, < 26(3 E,{T,^) , 
and Eq < El, so that we can rewrite the estimate of J2 as 

I J2I - J ||VA(^||i. < ^is'WVTWl, +CEliT, , 

with 

C := max (1, 26(3) 29/3 {2e^ + za) tfii'^il + tfii^) . 

Recalhng the estimate of Ji obtained at the beginning, in which we can 
bound ||V^||^2 by 2i^Ei{T,ip), we find that 

d^||Av,||i.+ e2||VA(/.||J. <c.e^El(T,^) + /ie2||VT||i.+Ci^;ll^,(^)^ 
where we have used the other shorthand 



II 1 2 

00 1 
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Finally, noting that ||T||x,2 < ||T||2,2 + ||T||/^2, and similarly for the derivatives 
of T, we obtain 

(52) d ^11 A(/.||i, + II VAv^ll^, < //el VT||i. + /i^l Vf ||i. 

Multiplying this inequality by z and adding it with (1471) . we can absorb 
||VT||^2 on the left-hand side provided that z^e^ is small enough compared 
to 59. Therefore, we eventually obtain ( HHl) by substituting a numerical 
constant times ^95e^ for z in the definition of C, here above. 

The next step consists in expanding -F(v^, Ay^, T)T, the first term in 
the right-hand side of (HSj) . and estimating each term by a power of Ei(T, ip). 
This is made in the following. 

Lemma 4.5 (final a priori estimate) In addition to the notations intro- 
duced in Lemma (see p. \2^, let us define the constants A, B, C, D, 
by 

A:=A^ + (3a{l + /i^/^)H^^ + 

C :=C° + eafx'/^W^\\f\\l, + 9a\W^\' + 9ae^fx'/^\\f\\%,, 
D:=D^ + e'{(39y/'a'Sfi (1 + a^^Vl/?^)'/') 
+ e(3fL{a'59y/'{l + e'fx%a'^9'SY/^) 

Then the first n modes (T = T + T,ip) = [T^ = T" + T, </?"■) solutions 
to fHT|l - fj42|) satisfy the energy estimate 

(53) ^^EliT, ^) + \59\\VT\\l, + a\\dt^\\l. + \e'60\V^^\\l. 

< AEliT, (/.) + D{B + El(T, ^))' + C, 

for all t G (0,t„). 

Proof: Recalling the definition of F in ( l33l) on p. [181 we have 
|F(y.,Av.,T)| < 2a|iy^|^ + a7z/;,iy^|r| + 2ae^ | Av;|^ + aeV;^ ^ ^ 

/l /2 /3 
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1. Since 



2\W'J^T < \W'J^ + W'^t"^ and 2\T\T < 3T^ + 



we readily have 

Jn 
hence 

(54) § [ hT<a^El{T,^) + c^, 

with 

2. By the Cauchy-Schwarz inequahty, we have 

/ |Ay,nT|<||A^||i,||T|U.. 
Jq 

Using the interpolation inequality (due to Holder) 

1/4|| ||3/4 



oo I 



\v\\t4 < \\v\\ A \\V 



Li ^ \\U\\l2 WVW^e , 



and the Sobolev inequality 



\h<\\v\\l2 + \\Vv\\l2 



(equivalent to the Sobolev embedding ^ already mentioned 
before), we thus infer that 

(55) / lA^plTl < ||A^||i.||T|U. + \\A^\\];,'\\TU4'^A^\\f. 
Jn 

As a consequence of and Young's inequality (Eq. (EI]) with a factor 
A to be determined afterwards, and again p = 4/3), we get 

9ljJ^<eae'\\A^\\l4T\\,.+9ae' (a^ || A^||i.||T||i. + A || VA^||i. 
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So, using that 

(56) \\A^\\l, < e'Safi9E,{T,^) , ||T||i, < peE,{T,^) , 
we arrive at 

(57) 9 [ hT<ahH^im"'E,{T,^f' 

Jn 

+ \^ah^Sfx(3^e^E^{T, iff + XOae^WVA^Wl^. 

We can now choose A for the last term in (1571) to be absorbed by 
VAv9||^2, the last term on the left-hand side of (jlH]) on p. l24l 
it suffices to take A small enough compared to aSfi. Hiding the multi- 
plicative constant in the < sign, we thus get from ( 1571) that 



(58) 
with 



e [ f,T-y'6fie\\VA^\\l,<d,El{T,^Y 



d, := e\f3ey/'a^fi (1 + aH' fi\f3df/^) . 



3. The way of estimating f^T is very similar to the one for /2T. 
Indeed, we have 

2 / |A^||T|T < / \A^\\W+ [ \A^\\f\\ 
Jn Jn Jn 

with 

/ |A^||f|^ < \\A^U4f\\l, 
Jn 

by Cauchy-Schwarz and the Sobolev embedding ^ L^, and 

/ lA^im^ <iiA^iu.|iTiii. + \\A^\\Lm]!"m\\f, 

Jn 

by exchanging the role of A(p and T in ( l55l) . Therefore, using again 
( 15T|) we have 

/ |A^||T|T < l|f + llAv^lli, + WA^hmh 
Jn 

+ X^A^\\U\T\\l, + X\\VT\\l. 
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for all positive A. Using again fl56p we thus infer that 

I f3T<9ae'^i'/^f\\%,+6a'^i'/'E^{T,^) 
Jn 

(59) + A^aeV^^llVTlli^. 

Now, if we want to let the last term in fl59p to be absorbed by the left- 
hand side term VT||^2 in f l48p . we choose A small enough compared 
to 6ae'^fi^/^. This yields 

[ fsT- i^^llVTlli. < 9ae'fi'/'\\f\\%,+6a'fi'/'E,iT,^) 
Jn 

(60) +ePfia^^\6ey/^Ei{T,ipf^ 

+ e^Pa^fi^eH^Ei{T,iff, 

hence 

(61) 9 [ f,T - l6e\\VT\\l,<c, + a3E,{T,^) + d3El(T,^f, 

Jn 

with 

and 

Finally, adding together the estimates in fl54p -f l58p -f l6ip and the energy esti- 
mate in fHHj) gives f l53|) with 

^ := + ai + ag , C := C° + ci + C3 , := + 4 + c?3 • 

Recall that T and ip in (153|) are actually and and of course depend 
on the rank n in the Galerkin approximation. However, their initial data 
To e L/^ and yjo ^ -f^n independent of n, and so is the initial energy 
-E'i(To, v^o)- By the energy estimate ( l53l) . there is a uniform time t* > of 
existence of T" and v?" - as solutions of the ODEs (l43l) - such that i?(T , 99") is 
bounded in L~(0,r), while ||VT 11^2, ||at(^''||L2 and ||VA(^"|| 1,2 are bounded 
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in L^(0,t*). Therefore, by the usual compactness theorems, up to extracting 
subsequences we have 

in L°^{0,t*;Hl) weak-*, in L'^{0,t*;H^) weak, in L'^{0,t*;L^) strong, and 
also almost everywhere, 

in L°^(0,r;L2) weak-*, _ _ 

in L°°(0,t*;L2) weak-*, in L'^{0,t*; H^) weak, in L^{0,t*;L^) strong, and 
almost everywhere. Furthermore, since 

dtif e L\0,t*;L'^) and if e L"^ {0 , t* ; H^) , 

we have by the Aubin- Simon lemma that 

G ^(0,r;i7^). 

Moreover, the definition of T" in ( 142|) shows that the sequence {dtT"') is 
bounded in the dual of L^(0,t*; ifg), which implies that 

T e ^(0,r;L2). 

Passing to the limit in (HTjl and ( 142|) we conclude that (T, ip) is a weak solution 
to (I32D on [0,r). 

Remark 4 (time of existence) 5y the energy estimate in (l53l) anc? an e/- 
ementary ODE argument we see that the time of existence t* is bounded 
from below by tl := l/(Z)£^i(To, v^o)^) (up to a multiplicative constant only 
depending on Q and the dimension d). This yields two comments. 

1. For fixed parameters a, P,'y,6,e,6 , by definition of Ei, t\ tends to in- 
finity when the initial data tend to a constant state {Tq = 0, (^q = 0) or 
(To = 0,ipo = 1). In other words, tl is all the more larger as the initial 
state is close to a single phase. 

2. We can also examine how tl varies with respect to the parameters a, (3, 
7, 5, e, 9 for fixed initial data. By inspection of all coefficients in the 
definition of D (p. [F^ we see that D is bounded on compact subsets of 

{(a, /3, 7, (5, £, ^) ; a > 0, /3 > 0, 7 > 0, 5 > 0, £ > 0, ^ > 0} , 
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and tends to zero in either one of the limits 

/3— 7-0 or ^ ^ or 6 ^ or 6 ^ oo . 

Note that, by definition of and 6 (see Eq. fl2U]) on p. 1/St = 
1/Pe = [55, so that St tends to infinity in either one of the first three 
limits, whereas Pe can tend to infinity (if P and 6 fixed), tend to 
zero (if 6 CO and (3 fixed), or can be kept bounded and bounded away 
from zero (if f3 — )■ and 5 — )■ oo with (3 and 1/6 of the same order). 
Unfortunately, it turns out that neither one of these limits is compatible 
with keeping i?i(To,</?o) bounded (which would imply t* — )■ oo), or at 
least D£'i(To, (/?o)^ bounded (which would give a uniform lower bound 
fort*), independently o/||To||l2 and ||A(/?o||l2. Indeed, we have 

2/3^^1 (To, y^o) > llTolli^ , 2«7>^|2^Ei(To,</^o) > e'6\\Aipo\\l., 
and looking closer at D we have 

D 0(/3i/2) , D 0(7^) , D 0(6'/') , D 0(9'/^) . 

So, the limit /3 — )■ penalizes ||To||l2, while 9^0 penalizes both 
||To||l2 and \\A(Pq\\i2, and both 7 — )■ and 5 — )■ 00 penalize ||A<y9o||L2. 
Of course, if we restrict our initial data to Tq = 0, the limit /3 — )■ zs 
allowed and thus gives t* — t- 00. This is no surprise because in this case 
we are basically left with the Allen- Cahn equation (the first equation in 
f in?]) with T = T), for which global existence is well-known. 

4.3 Continuous dependence on the data 

This final section is devoted to the proof of Theorem I4.3[ which gives both 
uniqueness and continuous dependence of the initial data. The tools are 
basically the same as in the existence proof (Theorem 14. 2p . namely, energy 
estimates and various inequalities of Sobolev and/or interpolation type. How- 
ever, the details are longer and somehow more technical. To gain some sim- 
plicity in the exposure we set all six parameters a, f3, 7, 6, e, and 9 equal to 
1 in (132 p . and the L°° norms of derivatives of W and u will be systematically 
hidden in the < sign. 

Step 1. We assume that (Tj, ipi), i G {1, 2} are two weak solutions to ( 132|) 
on [0,t*) for some positive t*. We are going to use repeatedly the notation 

[G] = G2-G, 
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for the difference of any two quantities G2 and Gi. The differences [T] and 
[if] then satisfy the system 

dtM-/\[^] + [W'{^)] = [u\^)Ti 
(62) { dAT] + [u\^)dM-m = [Fl 



By analogy with Eq defined in fl39|) we set 

eo(m,b])= / i[r]2 + i|v[(p]p + b][w^'(^)]. 



Similarly to the energy estimate (jH]) for Eq (except that the potential energy 
has been discarded here), we derive the identity 

(63) ^eo([T],M)+ [ \V[T]\'+ I \dA^f 

[u\v)T][dM-W{v)dM[T]+ I [dM[W'm + f [F][T]. 

Jn Jn 

Using the formula of differentiation 

[GH] = {G)[H] + [G]{H), (G):=^i±^, 
and the mean value theorem (for u') we obtain 

Wi^)T][dM - Wi^)dM[T] =[u'{^)m^)[T] - {T)[dM) 

(64) <l(5*^)ll[v^]|imi + |(T)||[v.]||[9,^]|. 

In order to estimate the first term of the right-hand side in (1641) . we can use 
the following inequality, whose proof is postponed to the appendix. 



(65) 



/ abc < e\\Vb\\l. + e\\Vc\\l. + C,(l + \\a\\l.){\\b\\l^ + \\4h). 
Jo, 



satisfied for all a G L^(f2), b,c ^ H^{Q), e > 0. Here above and in what 
follows, denotes a constant depending only on e, Q and d. The inequal- 
ity (^5^ gives 

l(5*V')llb]|imi<^||V[T]||i.+5||Vb]||i. 

(66) +c.(ii(5*^)rL. + i)(iimiii. + iivMiii.)- 
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To obtain a bound on the second term in the right-hand side of f l64p . we are 
going to use Agmon's inequahty ( ISUj) under the form 

(67) \\b\\U<\\b\\h + \\Vb\\,4Ab\\,., 
which imphes (by Cauchy-Schwarz) that 

(68) / \abc\ < e\\c\\l. + C,||a||i.(||6||i. + ||V6||i. + ||A6||i.). 
Jn 

We thus find that 

\{T)\M\\m\ < ewmwi. 

(69) + a\\{T)\U\M\l. + \\VM\h + IIAMIliO- 
Together with (166|) . (169!) gives 

(70) j[u\ip)T][dM - W{v)dM[T] < ^11 V[T] ||i. + e|| 

+ C,(1 + rn^nl. + ||(T)|i.)(l|[T]||i. + IIMIIi. + llVMIIi. + IIAMIliO- 

Let us introduce the shorthand Zq := 1 + \\ {dt(p) \\'^2 + ||(7')||^2, and 

ei([T], [^]) = eo([T], M) + f,\\A[cp]\\l. , e,{[Tl M) = ei([r], [^]) + M\\l,, 

where /i is a parameter that will be chosen later. The final estimate we are 
aiming at reads 

(71) ^e2<ze2, 62 = e2{[T],[ip]), 

for some 1 G L^{0, t*) depending continuously on the norms of Tj in L°°(0, t*ji^^) 
L2(0,r;L2) and the norms of in L°°(0, T; i/^)^ L\0,t*;H^), i G {1,2}Q. 
Once we have (17T]1 . the continuous dependence on the data as expressed 
in ( 15S]) will be clear because, by straightforward integration 



62 



(t) < e^o-W'^-e2(0), 



""^In particular, using the equation for (p in ()32|) . z may depend on the norm of dtip in 
L°°(0,r;i2) 
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and by definition 62 controls ||[<y9|||//2 and ||[T]||j;^2. 

The first step in the derivation of fl7T]) follows from the identity in 
and the inequality in ( !70|) . which together give 



(72) -eo([T], M) + ||V[r]||i. + \\dM\h < eWVmWl. + e\\[dM\\h 



+ C,,,Zoe,{[T],[^])+ / [dM[W'{v)]+ / [F][T], 

Jn Jn 

where zq e L^{0,t*), and C^,^ := max(l, 1/yu) C^. We will prove ([71]) by 
using an additional energy estimate for ||A[(^]||^2 and some estimates for the 
remaining terms 

'[dMlW'iv)] and / [F][T] 



in the right hand-side of fl72|) . 
Step 2. By gS]), 



(73) /[5t^][w^'(^)]<^ll[5*^]|li.+^illMlli- 

Step 3. There are several factors in the expansion of the term Jj^[F][T] 
that we bound successively and, to some extent, by decreasing order of dif- 
ficulty. 

By using the expansion formula 

[GHI] = (G) {H)[I] + {G) {I)[H] + {HI)[G], 

we obtain 
(74) 

[A^u'i^)T][T] = (A^) Wi^Wr + (A^) (T)K(v.)][T] + {u'i^)T)[Av][T]. 



By the inequality in fl65|l we have 

(75) / (A^) Wi^Wr < emml. + C,(l + ||(Av.)||t.)||[T]||i.. 
Jn 

and similarly 

jy{ip)T)[Av][T] < ^||V[T]||i. +e||A[y,]||i. 

(76) +C.(l + ||(T)||iO(imili. + ||Mlli.). 
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To estimate the remaining term in f l74p . we apply the inequahty (proved in 
the appendix) 

\abcd\ <e\\Vd\\l2+e\\Ac\\l2 

(77) + + llallt. + \\b\\l. + ||&||i.||V6||i.)(l|Vc||i. + ||rf||i.) 

to a = (A(/?), b = (T), c = [(p], d = [T]. This gives 



l(A^) (T)[z.'(^)][T]|<.||V[T]||i.+.||A[v.]||i. 

n 



+ + ii(A^)iit. + ii(T)iit. + ii(T)iii.iiv(T)iii.)(iivb]iii. + mr 



With fTfSl) and fl76l). we arrive at 



(78) / \[AW{m[T]\<e\\V[T]\\l2+MM\\h+ei). 

Jn 

Replacing A(p by W'{ip) in (JTHl) and using the estimate |[W^'(<yc)]| < l^c^lMI, 
we will obtain a similar estimate for the quadrilinear term |[lV((y9)z/'((/))T][T]. 
For the trilinear term 

[(A^ - W'{v)nT] = [\A^nT] + [W'iMT] - 2[A^W\^)][n 



by proceeding as in Step 2, we also have an estimate as in ( ITHl) . We thus 
conclude with Step 2 that 

(79) / I |[i^][T]|<e||V[T]||i.+^3(||b]||i.+ei). 
Inserting (179|) in (!72|) gives 



(80) ^eo + ||V[T]||i. + ||a4v.]||i. < e||V[T]||i.+.||[a,v.]||i.+^4(||Mlli.+ei). 

To eventually obtain flTTj) . we still have to estimate — 1| A[y9]||^2. 

Step 4. Taking — A^[(y9] as a test functioij^ in the equation satisfied by 
[ip] in fl62l) . we obtain 
(81) 

^||A[v.]||i. + ||VA[v.]||i. < ^ |[V(z.'(^)r)][VA^]| + jj[VW\v)][VAv]\. 

^To be exact, we take — A^[(^]" as a test function, where [tp]" := X^iLid"^]' V'O'Pi is the 
Galerkin approximation of rank n to and, in a final step, pass to the hmit n — > +oo in 
the estimates obtained for [^p]^ . 
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We have 

(82) / |[Viy'M][VA(/,]|<5||[VAv.]||i.+C,||[^]||i.. 
Jn 

Besides, we have the expansion 

(83) I [V{u'{^)T)] = {u'{^))[VT] + (VT) + {Vu'{^))[T] + (T) [Vu'{^)] 
which has to be multiphed by [VA</)]. We have 

(84) / \{u'{^))[VT][VA^]\<s\\[VA^]\\l, + Cm\l^ 
Jn 

for the first term, and by (!68|) 

\{VT)[u\^)][VA^]\<e\\[VA^]\\l, 



^5) +Ce||(Vr)||i.(||[^]||i. + ||VMIIi. + ||A--'2 



L2y 



for the second term. To estimate the third term, we apply the inequahty 
(proved in the appendix) 



'a\\l. + \\a\\l.)\\b\\l. 



(86) / \abc\ < s\\c\\l2 + Ce\\Vb\\l. + Q(l + || Vc 
Jn 

to a = (Vv?), b=[T], c= [VAip]. We obtain 

{\/u'i^))[T][VA^]\ < e||[VA(^]||i. + C,||V[T]||i. 



^7) +C,(l + (||V(^||i.r+(||Av.||i.))|[T]||i.. 



Since [Vi^'((/?)] = {iy"{ip))[Vip] + {'Vip)[iy'{(p)], the fourth term can be spht 
again. Similarly as in f l85|) . we have 



\{T){u"{^))[Vcp][VA^]\<e\\[VA^]\\l, 
(88) + Q||(T)||i.(||M||i. + iivMiii. + ||A[v.]||i.; 

There remains to bound 

'l(T) (V^)M[VAv.]|. 
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To this purpose, we apply the inequahty (proved in the appendix) 

/ \abcd\ < e\\d\\l, 
Jn 

(89) + Q(||a|U. + II Va|UO(ll&lli^ + imihMcWl^ + || Vc||i. + || Ac||i.) 

to a = (T), 6 = (Vifi), c=[(f],d= [VAcp]. We then deduce from (|8T|), (|82|), 
(IMj) . dSSD, dHZD, dHSD, taking e small enough, the estimate 

(90) ^IIAMIIi. + lljVAMIIi. < C4VT\\l, + CzEimWh + ^i)- 

Choosing /i so that fiCi is small enough, it then follows from ( IHOl) and (l90i) 
that 

(91) ^ei + ||V[T]||i. + ||9,[(/.]||i. < 5||V[T]||i.+5||[a,(^]||i. + ^6(||M|li.+ei). 

To conclude, we perform an energy estimate on [if]. Multiplying the first 
equation in (l62l) by [ip] we obtain 

^IIMIIi^ + livb]||i. < c (^MWl, + ||[r]||i.) + JjiTMU) 

<c(IMIIi^ + limiliO + IK3^)llLHIMIII) 

< c(i + ll(T)||iO(IIMIIi. + limili.) + llvMlli.)- 

Adding with ( 19T]) we get 

^(ei + l|[^]|liO <^r(ei + ||M||i.), 
which is the claimed estimate ( 1711) . This completes the proof of Theorem 14.31 

Appendix 

We give here the proof of various functional inequalities used previously. 
Proof: [Proof of ([65])] 

/ abc < e\ml. + e|| Vclli. + C,(l + ||a||i.)(ll&lli. + l|c||i.) 
Jn 
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for n C W'-, d < 3. By Holder's inequality, 

abc < \\a\\L2 \\b\\L6 ||c||l3. 



In addition, we have the interpolation inequality ||c||j;^3 < ||c||]^(^||c||]^{^, and 
applying twice Young's inequality (15T]) (once with p = 4/3, once with p = 
3/2), we get 

ll^'lUa • llcll^/^ ■ ||a|U.||c|lK' < e\\b\\l, + e||c||ie + C,||a||i.||c||i.. 

We can conclude thanks to the Sobolev inequality ||&|||6 ^ 11^11^2 + ||V6||^2- 
Proof: [Proof of ([77])] 



\abcd\ <e\\Vd\\l2 + e\\Ac\\l2 

+ + MI2 + ml. + ||6||i2||V6||i.)(l|Vc||i. + \\d\\l.] 
for Q cMf^, d < 3. By Holder's inequality, we have 

/ \abcd\ < ||a||L2||fe||i3||c||Loo||(i||i6, 
Jn 

where we can apply Agmon's inequality ( 1671) to c: 

llclUoo <C(||c|U2 + ||Vc||i/'||Ac||^/'), 
and the interpolation inequality II&IIl^ < ll^llii^ll^llle^- This gives 

\abcd\ < C{\\ch. + ||Vc||i/^||Ac||^/^)||a|U2||6||^/^||fe||^/^|M|Ua. 
By the Sobolev inequality \\b\\L(i < C(||6||l2 + ||V6||l2), it follows that 
\abcd\<Ci\\c\\L2 + \\Vc\\'/,'\\Ac\\]^,^) 



(92) X \\ah4b\\'/'{\\bC' + ||V6||^/^)(||d|U2 + ||Vrf|U2). 

The term with the highest order derivatives in f l92|) is 

||Vc||^/^||Ac||^/^||a|U2||6||i/^||V6||l/^||Vrf|U2 = 

IIV^L- ||Ac||^//- ||Vc||^/^||a|U2||6||i/^||V6||i// 
<.||Vrf||i.+5||Ac||i.+C,||Vc||i.||a||i.||6||i.||V6||i2 
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by Young's inequality in the form already used above, namely 

We finally obtain fl77|) by using similar estimates for the lower order terms 
in ([92]). 

Proof: [Proof of §6^] 



\abc\ < e\\c\\l, + C,\\Vb\\l, + + || Va||i. + ||a||i.)ll&lli. 

'n 

for Q C M.^, d < 3. As above, applying successively the Holder, Young, 
interpolation, and Sobolev inequalities give 



/ \abc\ < \\c\\L2\\b\\L3\\a\\L6 
Jn 



< e\\c\\l, + CMlMl^ 

<E\\c\\l, + CML4b\\L44l<^ 

< 4c\\h + CMI^ + CsMUMh 

< e\\c\\h + C^VbWl. + + llValli. + ||a||i.)ll&lli- 
Proof: [Proof of (Ea])] 

\abcd\ < e\\d\\l2 

+ CMah^ + II ValUOdl^lli^ + II V6||i.)(l|c||i. + || Vc||i. + || Ac||i.) 

for Q C M'', d < 3. Again, by the Holder, Young, interpolation, Sobolev, and 
Agmon inequalities, 

\abcd\ < \\d\\L2\\a\\L3\\b\\Le\\c\\Loo 
<dM||?, + a||a|'2 11^.112 



1^2 -T '-^£||W||L3|niL6 1l'-|lL°° 

< e|M||i.+C,(||a|U. + llValUOdl&lli^ + llVfelli^) 
X (||c||i2 + ||Vc||i2 + ||Ac||iO- 



Note added in proof 

After completing the manuscript, we became aware of a recent work by 
Feireisl, Petzeltova and Rocca [7], in which they analyze a closely related 



40 



phase transition model with microscopic movements. In particular, they 
take into account a quadratic term |(9t(y9p [ip being the order parameter) in 
the equation for the temperature, like the second equation in our system 
f lT6|) . which quite equivalently involves the quadratic term A(pdtip. They do 
not address the sharp interface limit, but they study local existence, and 
even global existence of a certain type of weak solutions, obtained by passing 
to the limit in a regularized system, in which the diffusion term —AT is re- 
placed by — div ((1 + £:T'")VT), m large enough. Global existence of classical 
solutions, and even of weak solutions according to the notion we have used 
in the present paper (see Definition 14. ip remains an open question up to our 
knowledge. 
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